Spin-orbit coupling in ferromagnetic Nickel 
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We use the Gutzwiller variational theory to investigate the electronic and the magnetic properties 
of fcc-Nickel. Our particular focus is on the effects of the spin-orbit coupling. Unlike standard 
relativistic band-structure theories, we reproduce the experimental magnetic moment direction and 
we explain the change of the Fermi-surface topology that occurs when the magnetic moment direction 
is rotated by an external magnetic field. The Fermi surface in our calculation deviates from early 
de-Haas-van- Alphen (dHvA) results. We attribute these discrepancies to an incorrect interpretation 
of the raw dHvA data. 
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The limitations of density functional theory (DFT) 
when treating the electronic and magnetic properties of 
transition metals become evident most clearly in the 
case of Nickel. The DFT cannot reproduce gross fea- 
tures such as the width of the 3d-bands (4.5 eV versus 
3.3 eV experimentally^ 2 ^) , nor important details such 
as the exchange splitting. The exchange splitting in the 
DFT is almost 0.7 eV and rather isotropic over the Fermi- 
surface, whereas, experimentally, it is found to be much 
smaller and strongly orbital dependent: A e ss 0.17 eV 
and At 2g rs 0.33 eV. As a result, even the Fermi surface 
topologies do not match, because of the position of the 
X%,[ energy: above Ep in DFT, yet below E-p experi- 
mentally; thus only one hole ellipsoid exists around the 
X point, versus two in DFTi^£ 

Even more limitations of DFT become evident when 
the effects of the spin-orbit coupling are considered. The 
magnetic anisotropy energy has the wrong sign for Nickel 
(and for Cobalt), while it has the correct sign for Iron, 
yet is too small by a factor of three£ In Nickel, the easy 
axis is along [111] and approximately 3/^eV per atom 
are needed to rotate the magnetic moment axis into the 
[001] directional^ Moreover, a detailed low-temperature 
study of the magnetic anisotropy constants K\ , Ki , K$ 
by Gersdorf^ has revealed a change in the Fermi-surface 
topology when the magnetic-moment axis is rotated into 
the [001] direction: A small second hole ellipsoid appears 
around the X(001) point, but not around the X(100) and 
X(010) points, now inequivalent to AT(001), because of 
the underlying tetragonal symmetry. 

It is generally accepted that the discrepancies between 
the DFT and the experimental results are mainly caused 
by an insufficient treatment of the electronic correlation 
in an effective one- particle theory. In the past, all at- 
tempts to combine the DFT with more sophisticated cor- 
related electron theories have only led to partial improve- 
ments of the results for Nickel; see, e.g., the GW approx- 
imation in Ref. |{j. 

In a recent work^ we were able to show that a general- 
ized Gutzwiller theory provides a consistent picture of the 
quasi-particle band-structure of Nickel. Neglecting spin- 
orbit coupling, all basic problems of the DFT calculations 
on Nickel have been resolved. Our theory employed ap- 



proximately 2 10 variational parameters representing the 
occupancies of all atomic multi-electron states within an 
open 3d shell (see below). 

In this letter we present results for the case when spin- 
orbit coupling is included. In order to cope with this com- 
plication, the Gutzwiller theory had to be extended^ to 
allow for rotations in the eigenvector space of the atomic 
atomic multi-electron states, resulting in many more vari- 
ational parameters. Employing this generalization we ob- 
tain the correct magnetic anisotropy energy, and, more 
importantly, reproduce the change in the Fermi-surface 
topology found by Gersdorf. 

To investigate transition metals we start from multi- 
band Hubbard models of the general form 



H = £ t 



a, a' a| 
i,j C hv C j 



](T / + 2^ H\oc,i — H o + H\ oc 



(1) 



Here, the first term describes the hopping of electrons 
between spin-orbital states a, a' on lattice sites i,j, re- 
spectively. The Hamiltonian 



m 



loc,', 



Hc.i + H c a + Hi 



SCM 



(2) 



contains all local terms, i.e., the two-particle Coulomb 
interaction Hc,i, the crystal field energies H c fi and the 
spin-orbit coupling -ffso.i- I n the case of Nickel, we work 
with a basis of 3d, 4s, and Ap orbitals. 

We have determined the bare hopping-parameters 
in the one-particle Hamiltonian Hq and the crystal- 
field energies in H c { by means of a tight-binding fit 
to the paramagnetic DFT band structure 4^ Due to 
the large band- width of the 4s and 4p bands, only the 
Coulomb-interaction within the 3<i-shell is taken into ac- 
count. The spherical approximation is used, i.e., we ex- 
press the Coulomb interaction through the three Racah- 
parameters A, B, and CJ^- Note that cubic site symmetry 
would allow ten independent interaction parameters. In 
order to reproduce the experimental d-band width in our 
approach, we need a Racah-parameter A ps 9eV. The 
Racah-paramctcrs B and C are assumed to be close to 
their atomic values^ 1 , B m 85 meV and C rj 400 meV, re- 
sulting in a value J of J = 7B/2 + 7C/5 « 0.85 eV. The 



spin-orbit coupling parameter £ in the spin-orbit Hamil- 
tonian 

Hso,i = /]n(<r\lx<rx + ly°y + lz^zW)ct iCT Ci^' (3) 

is chosen as C = 80meV. Note that the Hamiltonian © 
only contains d-orbitals. 

In the Gutzwiller theory, the following Ansatz for a 
variational wave- functio n 10 i 13 ' 14 is used to investigate the 
multi-band Hubbard model ((T|) 



|*G> = Aj|*0> = n^l*0>- 



(4) 



Here, |^o) 1S a normalized single-particle product state 
and the local Gutzwiller correlator is defined as 



A = X) A r>|r>«<r / 



(5) 



r,r' 



The states \T)i form some arbitrary basis of the atomic 

Hilbert-space and the (complex) numbers Ai? r , are varia- 
tional parameters. For Nickel, we work with a correlation 
operator |(SJ) in which the states |r)j are the eigenstates 
of the atomic Hamiltonian Hc,i- The non-diagonal ele- 
ments of the variational parameter matrix Ar,r' are as- 
sumed to be finite only for states |T), |r") which belong 
to the same atomic multiplet. This is consistent with 
the spherical approximation for the Coulomb-interaction. 
In the case of Nickel, it is sufficient to work with non- 
diagonal parameters Ar,r' in the d 7 , d a and d 9 shells. 

The expectation value of the Hamiltonian fTJ) can be 
calculated analytically for the Gutzwiller wave function 
(fj| in the limit of infinite spatial dimensions^ We use 
the exact results in this limit as an approximation for our 
three-dimensional model system. The Gutzwiller energy 
functional can also be obtained in slave-boson mean-field 
theor y 16 i 17 In infinite dimensions one finds 

(-Hioc,i}* G = 2_, A r2Ti^r 3 ,r4 £, r 2 c ,r 3 m ri,r4 ( 6 ) 
ri...r 4 

for the expectation value of the local Hamiltonian in |T]) , 
where 



^T° c ,r 3 = i^2\Hi OCt i\T 3 ) . 

< lt r< = ((|ri>(r 4 |))* 



(7) 
(8) 



The local expectation value (|SJ) is readily calculated by 
means of Wick's theorem. For the expectation value of a 
hopping operator in (|TJ) one finds 



\S,<Ji°J,(T2' 



)*G = J2 & \f^) ( B l<T[Cj,<T' 2 )*0 ■ ( 9 ) 



The renormalization matrix q a a in ^ can be calculated 
most easily when an orbital basis is used which has a 
diagonal local density-matrix with respect to \^o), 



a. 



- /zi 



<cLc i]O .,)* =5 a 



(10) 



If C° a , is non-diagonal for a one-particle product state 
|^o) one can always transform the orbital basis in order 
to ensure that Eq. (JTQJ) holds. In the case of a diagonal 
local density-matrix (|I0[) . the renormalization matrix in 
® reads 



£' = ^ £ A 



r 2 ,r,-V 



,r 4 (r 2 |4ir 3 > 



(11) 



ri...r 4 



x <(|r 1 )<r 4 |c CT 0)* o . 

where, again, the expectation value with respect to |\&o) 
is calculated with Wick's theorem. 

The variational ground-state energy must be mini- 
mized with respect to the variational parameters Ar,r' 
and the one-particle product wave- functions |^o)- H 
has been foun d 10 ' 18 that the optimum state |\&o) is the 
ground state of the effective one-particle Hamiltonian 






A.,'- (12) 



Here, we introduced the renormalized hopping matrix el- 
ements 



t: 



1,3 



E^MY 1 ®* 



(13) 



and the Lagrange-parameters r) a ,ir' which are used to op- 
timize the energy with respect to the local density ma- 
trix (fTU|) . Within a Landau Fermi-liquid approach one 
can further sho w 10 ' 18 that the eigenvalues E 1 (k) of Hq S 
are the quasi-particle excitation energies that can be 
compared, for example, to ARPES experiments. Most 
important for the quasi-particle band-structure are the 
Lagrange-parameters rj^ a , for the d-orbitals. The two 
(diagonal) Lagrange parameters for the s and p orbitals 
are adjusted in order to fix the total d-electron number J£ 
The inclusion of spin-orbit coupling in the Gutzwiller 
theory complicates the numerical minimization signifi- 
cantly. Both the d-part of the local density-matrix (fTU]) 
and of the hopping renormalization matrix (|If [) are no 
longer diagonal. The number n- m of independent ele- 
ments depends on the magnetic moment direction, we 
find n io = 22 for /Z||[lll] and n ic = 18 for #||[001]. As 
a consequence of the reduced symmetry, we could work 
with up to ni c independent d-shell Lagrange parame- 
ters vfc a , in order to minimize the total energy. Nu- 
merically, however, such a minimization would be quite 
costly since each variation of these parameters involves 
many momentum-space integrations. We therefore work 

with a simplified effective Hamiltonian Hq S that contains 
effective parameters only for all physically relevant one- 
particle terms. 

In cubic symmetry, there exist only four independent 
matrix elements of the local (d-electron) density-matrix. 
The trace of the matrix is fixed by the total d-electron 
number. The three remaining matrix elements are gov- 
erned by parameters 77^ a , which are given by the orbital- 
dependent exchange fields A t2g , A Cg and the effective 
crystal- field splitting e^p. 
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FIG. 1: Quasi-particle band structure along the A-line around 
the X z -point of the Brillouin zone for magnetic-moment di- 
rections jl || [001] and /2 || [001]. The inset shows an enlarged 
view of the band structure around the Fermi energy which dis- 
plays the additional hole ellipsoid for /2 || [001] more clearly. 



The non-cubic symmetry resulting from the addition 
of the spin-orbit coupling adds many more formally inde- 
pendent ^cr' terms. Both for /u||[001] (tetragonal sym- 
metry) and for /x||[lll] (trigonal symmetry) there are 
two more exchange-fields and two more crystal-field split- 
tings. All these eight terms are included in our simplified 

effective Hamiltonian Hq S . In the spirit of the spherical 
approximation 1 ^, a Hamiltonian H§q is included in Hq S 
that has the same form as Hso only with £ replaced by 
£ off . As a result, we have to minimize the total energy 
with respect to nine 'external' parameters in our simpli- 
fied Hamiltonian Hq S ■ 

The numerical minimization is much more time-con- 
suming for a system with spin-orbit coupling than with- 
out. First, spin-orbit coupling requires the momentum- 
space integration to be extended from l/48th to the full 
Brillouin zone. Furthermore, the small values of the 
anisotropy energy necessitate a much finer mesh for the 
momentum-space integration. Second, the energy needs 
to be minimized with respect to nine external variational 

parameters in Hq S . Altogether, the minimization of the 
total energy is approximately 10 4 times more time con- 
suming for a system with spin-orbit coupling than in the 
absence of Hso ■ 

We carried out the minimization of the variational en- 
ergy with respect to the 'internal' parameters Ar,r' an d 
the external parameters for both magnetic moment direc- 
tions /x| | [111] and /2||[001]. The optimum value of the ef- 
fective spin-orbit coupling is Ccff ~ 68 meV in both cases, 
about 15% smaller than the bare value £ = 80meV. 



FIG. 2: Same as in Fig.[T]but around the X^-point. Note that 
the X2|-band is below the Fermi energy for both magnetic- 
moment directions. 



There seems to be no simple rule that determines the 
relative size of £ and £ e ff . For example, for Iron we found 
an effective spin-orbit considerably larger than the corre- 
sponding bare value. In our calculations for Nickel, the 
anisotropy energy is -E aniso « 3.5 /xeV per atom, quite 
close to the experimental value -E exp ~ 3.0 ^eV. Note 
that this energy difference has to be calculated quite care- 
fully within the Gutzwiller approach. In particular, one 
has to keep in mind that any approximation on the pa- 
rameters Ar,r' that reduces the variational flexibility may 
lead to a grossly overestimated anisotropy energy. This 
is a serious problem, in particular, in the case of Iron. 
For Nickel, however, the active multiplet states belong 
mostly to d 8 and d 9 . Here, a mixing of states |T), |T') 
has little effect on the variational energy, and even a di- 
agonal variational parameter matrix Ar,r' ~ <^r,r' would 
lead to reasonable results. 

In Figures Q] and [5] we show the quasi-particle band- 
structure that arises from our calculation around the Ap- 
points X z = (001) and X x = (100). When the magnetic 
moment is along the easy axis, the band-structure around 
both Appoints coincides and the minority state X2J. is 
below the Fermi-energy^ For a magnetic moment along 
the [001]-direction, however, the two states A^j. have dif- 
ferent energies. The X21 state at X x remains below the 
Fermi level, whereas the corresponding state at X z cre- 
ates a new hole pocket around this AT-point. This is the 
scenario proposed by Gersdorf. 8 

In Figure [3] we show Fermi-surface cuts that we find 
within our Gutzwiller theory. The experimental values 
are taken from dHvA experiments by Tsu£ and by Stark 
as reported in Ref. [l5|. The agreement is quite satisfac- 
tory along high-symmetry lines, whereas there are signif- 
icant discrepancies away from them. We do believe that 



u x z X x u 




FIG. 3: Fermi-surface cuts with various planes in the Brillouin zone. Lines: Gutzwiller theory including spin-orbit coupling; 
Squares and triangles: experimental data reported infla: Dots: experimental data of Tsu£. 



the wiggles that appear in the experimental data are, in 
fact, spurious. The derivation of a Fermi surface from 
the raw dHvA data requires an expansion in Fermi sur- 
face harmonics, with the coefficients of the harmonics to 
be determined form least squares' fits to the data. Possi- 
bly, an over-determination occurred which led to unphys- 
ically large higher harmonics coefficients and resulted in 
the wiggles. We propose to redo these measurements. 

In summary, we have resolved the long-standing prob- 
lem to explain theoretically the electronic and magnetic 
properties of elementary fcc-Nickel. Our calculations are 
based on the Gutzwiller variational theory which is a 
powerful tool for the investigation of Fermi-liquid sys- 



tems with medium to strong Coulomb interaction. For 
such systems, state of the art band-structure theories 
usually fail. Our results for the quasi-particle bands are 
in very good agreement with ARPES experiments and we 
find the experimental Fermi-surface topology. Further- 
more, we are able to explain the subtle effects that the 
spin-orbit coupling has in Nickel. Our theory yields the 
correct anisotropy energy and we confirm the Gersdorf 
scenario: The Fermi-surface topology changes around the 
X-point (001) when the magnetic moment direction is 
rotated from /x||[lll] to /2||[001] by an external magnetic 
field. 
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